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1.0 Experiment Site and Schedule 

The Colorado State University surface observing systems described in this document 
were deployed at the Parsons’ Kansas site dining the FIRE II Cirrus Special 
Observing Period (SOP) from 13 November to 07 December 1991. The specific 
location of most of the CSU instrumentation was at a recreational area maintained 
by Kansas Gas and Electric Co.; this site was approximately eight miles east and 
two miles south of the town of Parsons. The geographical coordinates of this site 
are 37*18 min. N. latitude and 95* 07 min. W. longitude; site elevation was 269 
meters. In addition one surface meteorological and broadband flux observing site 
was maintained at the Tri City Airport which is approximately eighteen miles due 
west of Parsons (37 * 20 min. N. latitude, 95 * 30 min. 30 sec. W. longitude. Figure 
1 presents a map of the locations of the CSU deployment sites. 

At the main Parsons’ site the instrumentation was located directly adjacent to and 
north of a recreational lake. Under most cirrus observing conditions, when the wind 
had a significant southerly component, the lake was upwind of the observing site. 

Table 1 lists the measurements and observations collected during the experiment. 
These measurements may be grouped into five categories: Surface meteorology; 
infrared spectral and broadband measurements; solar spectral and broadband 
measurements; upper air measurements; and cloud measurements. 

Table 2 presents a summary of observations collected at the Parsons site during the 
SOP. One should be aware that the wind profiler, the laser ceilometer, surface 
meteorology and surface broadband radiation instrumentation were operated on a 
continuous basis. All other systems were operated on an "on demand" basis when 
cloud conditions merited the collection of data. 
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Figure 1: Map of the locations of CSU deployment sites. 




Table 1 Instrumentation deployed at FIRE II Cirrus Parsons site. 


Category 

Instrumentation 

Observed Variable 

Surface Meteorology 

Campbell Surface 
Station 

T, T d , u, v 

Infrared Radiation 

Pyrgeometer 

3-50 /xm flux 


PRT-6 

11 /xm radiance 


Interferometer 

5-15/xm radiance 

Solar Radiation 

Pyranometer 

.3-2.8 /xm flux 


Pyrheliometer 

.3-2.8 /xm direct component 
.53 jim long pass filter 
.695 /xm long pass filter 
1.0 /xm long pass filter 


MFOV 

.4-1 .0 /xm;5, 10,20, 28 * FOV 


Sun Photometer 

.38,.412,.5,.675,.862 /xm 

Upper Air 

Radiotheodolite 

P, T(p), T d (p), u(p), v(p) 


403 MHz Wind 
Profiler 

RASS 

u(z), v(z), w(z) 
T(z) 

Cloud Measurements 

Laser Ceilometer 

Video All Sky 
Camera 

Cloud base height 
(z< 7.58km) 

Cloud cover, cloud type 


Cloud lidar 
(Operated by 
Joe Alvarez, NASA 
Langley) 

Cloud base, thickness, 
Backscatter coefficient 
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Surface Meteorology 

2.1 Campbell Surface Station 

The radiation and meteorological data collections employed data loggers 
from Campbell Scientific Inc. Each remote station was slightly different from 
the others in order to make use of the systems on hand. 

Station 1 was located on a dyke approximately 120 feet South of the wind 
profiler location at the Parson’s site. This system consisted of a Campbell 
CR10 data logger with an AM32 thirty-two-channel multiplexer to increase 
the number of effective input channels. The meteorological instruments 
consisted of an R.M. Young model 05103 wind monitor and a Campbell 
model 207 air temperature/relative humidity sensor. The radiation sensors 
included Eppley PSP pyranometers for hemispheric solar and near IR 
readings, an Eppley PIR pyrgeometer for hemis ph eric IR data, and two 
Eppley NIP pyrheliometers-one filtered, one unfiltered-for direct solar 
measurements. The hemispheric instruments were mounted looking upward 
from the h-frame support structure. The pyrheliometers were mounted on 
Licor solar trackers. 

Station 2 was located at the Parsons Tri-City Airport, approximately 25 miles 
West of the profiler site. A Campbell CR21X data logger powered this 
system. The meteorological and hemispheric radiation sensors were the same 
as on Station 1 . Measurements of the direct solar component were not made 
at this site. 

Station 3 was used to collect sun photometer data. Another Campbell CR21X 
data logger received data from CSU’s Multiple Field of View (MFOV) 
photometer, and from a Mainz Sunphotometer (model MZ II 585). Both 
photometers were also mounted on Licor solar trackers. 

All stations collected instantaneous data at two-minute intervals. Station 3’s 
data rate was occasionally increased for calibration and comparison purposes. 

Data were downloaded from stations 1 and 3 each morning, and from station 
2 every two days. The data were initially stored on 740 kbyte floppy disks 
in a comma-delineated ascii format. This data was archived to Exabyte tape 
from a 386 PC using the Novaback backup utility, and is also archived under 
the FIRE data account on the HP9000-730 workstation. 
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2.2 Temperature and Relative Humidity 



Figure 2: Surface Instrument: 

Lower left: Filtered and unfiltered pyrheliometer; Center: 

Model #207 temperature and relative humidity unit; Upper 
right: Wind speed and direction sensor. 

The Campbell Model 207 unit measures temperature and relative humidity 
using two sensors combined into one probe. Separate CR7 input/output 
instructions are used to read each one. The probe was mounted on an H 
frame with a radiation shield. 

The probe contains a Phys-Chemical Research PCRCII RH sensor and a 
Fenwal electronics UUT50J1 themistor configured for use with CSrs CR7 
data logger. 

The instruction used to read the temperature sensor provides a 4V AC 
excitation, makes a single ended measurement, and linearizes the result with 
a 5th order polynomial to output temperatures in * C. 
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The instruction used to read RH sensor provides a 3V AC excitation, makes 
a single ended measurements and linearizes the result with a fifth order 
polynomial to % RH. It performs the required temperature compensation 
before outputing the result. 

The following are the linearization errors: 

thermistor 


-40 to 56' C 

< 1*C 

-36 to 53 *C 

<.5*C 

-33 to 48 *C 

<.1*C 

sensor 

12 to 100% 

<3% 

25 to 94% 

<1% 


The temperature compensation for the relative humidity is given by: 
RH = RH + .36(25-T) 


where T - temperature of the air in deg C. 

The accuracy of the measurements for the temperature sensor depend on a 
combination of FenwaPs interchangeability specification, precision of bridge 
resistors and linearization errors. The overall accuracy is ±.2*C. 

The accuracy of the RH sensor is a combination of the RH chip, precision 
resistors, and linearization error. The overall accuracy is quoted as better 
than 5% in the range 12% to 100%. 

For additional information on the data gathering system, data storage or data 
location, see section 2.1. 

2.3 Wind Direction and Speed 

Wind Speed Measurements: 

Range 
Sensor 
Pitch 

Distance constant 
Threshold sensitivity 
Transducer output 

20 V p-p at 1200 rpm 


0 to 60 m/s (80 m/s gust survival) 
18 cm 4 blade helicoid propeller 
29.4 cm 

2.7 m for 63% recovery 
.9 m/s 

100 mV p-p at 60 rpm 
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Wind Direction Measurements: 
Range 
Sensor 

Damping ratio 
Delay distance 
Threshold sensitivity 

Transducer 


360* mechanical, 355* electrical 
Balanced vane 38 cm turning radius 
.25 

1.3 m for 50% recovery 
1 m/s at 10* displacement 
1.5 m/s at 5* displacement 
.25% linearity-rated 1 watt at 40 * C, 
0 watts at 125* C 


Wind speed and direction were monitored using a propeller anemometer 
manufactured by R. M. Young. The main housing, nose cone, propeller, tail, 
and internal parts are injection-molded, uv-stabrlized plastic. The propeller 
rotation produces an AC sine wave signal with frequency proportional to 
wind speed. Three complete sine wave cycles are produced for each 
propeller revolution. 

Vane position is transmitted by a 1 OK ohm precision conductive plastic 
potentiometer. The output signal is an analog voltage directly proportional 
to azimuth angle. 

The monitor was mounted on an H frame stand and is calibrated before 
shipment. It requires no adjustments. 


References: 

Campbell Scientific, 1983: Model 207 Temperature and Relative Humidity Probe 

Instruction Manual. 

21X Micrologger Instrumentation Manual. Revision: 21X.1M. Copyright, 1984. Campbell 
Scientific, Inc. 

CR 10 Measurement and Control Module Operations Manual. Revision: 11/88. Copyright 
1986, 1988. Campbell Scientific, Inc. 

R. M. Young Co., 1988: Wind and Temperature Instruments. 
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3.0 Infrared Radiation 


3.1 Pyrgeometer 



Figure 3: Surface irradiance sensors: 

Left: .695-2.8 micrometer pyranometer; Center: .3-2.8 

micrometer pyranometer; Far right: 3-50 micrometer 

pyrgeometer. 

3.11 Instrument Description 

The Eppley Precision Infrared Pyrgeometer measures irradiance in the 4- 
100 /im spectral region of terrestrial radiation. It uses a blackened multi 
junction thermopile to create a voltage which ideally is linearly related to the 
net gain in radiant power. The thermopile is coated with Parsons black 
lacquer and emits like a black body at a temperature T. A thermistor is used 
to measure the sink temperature of the thermopile (see figure below). 

This pyrgeometer achieves its aim by isolating the desired spectral region 
with a single dome of silicon used with an interference filter vacuum 
deposited on the inside of the dome. 


3.12 Specifications 

Sensitivity - 5 /iV/Watt/m 2 

Impedence - 700 ohms 
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± 2% (-20* C to 40* C) 
± 1% (0 to 700 watts/m 2 
2 seconds 


Temperature Dependence 
Linearity 
Response Time 

Cosine Response - < 5% from normalization 

(Insignificant for a diffuse source) 

The orientation of the instrument has no effect on performance, but the 
accuracy of the instrument is very dependent on calibration which requires 
special care due to the dome-sink temperature differences. 

3.13 Data Reduction 

Using the reduction equations from Albrecht and Cox, 1976, the incident 
radiance, L, upon the pyrgeometer is given by the following equation. 

L = K, E + oeJ s A - K2 ff(' T d 4 - T s 4 ) 

T g - sink temperature 

T d - dome temperature 

L - incident radiation 

E - voltage 

The constants K, and for the two stations are given below. 

Station 1: K, = -243.309 iq = 3.85 

Station 2: K, = -234.742 = 4.1 

For any additional information concerning data gathering systems, data 
storage or data location, see section 2.1. 

Data are archived at NCDS. 

References: 

Albrecht, B. A. and S. K. Cox, 1976: Pyrgeometer data reduction and calibration 
procedures. Department of Atmospheric Science, CSU, Paper No. 251, 1-16. 

The Eppley Laboratory: An instrumentation for the measurement of the 

components of solar and terrestrial radiation. 


11 



3.2 PRT-6 Infrared Radiometer 


» 



Figure 4: PRT-6 infrared bolometer (11 micrometers). 


3.21 The Instrument and Its Purpose 

The PRT-6 is an all-purpose chopped bolometer which can passively 
sense ambient infrared targets within the range of 2 to 20 microns. 
It can be deployed vertically or at any zenith angle orientation, and 
at a field of view of either 2 or 20 degrees. The output is a voltage 
signal, recorded every 5 seconds, which possesses a linear relationship 
to incident irradiance. 

3.22 Capabilities and Specifications 

Throughout the entire 1991 FIRE IFO Experimentlhe instrument was 
set at a field of view of 2 degrees and made use of a filter whose 
transmission curve is displayed in Figure 5. This filter effectively 
narrowed the spectral band to a range of about 885 to 945 inverse 
centimeters. Most of the measurements were made using a vertical 
orientation, although zenith angles of 15, 30, 45, 60, and 75 degrees 
were also utilized. 


During the experiment an electrical attenuation component was 
sometimes added to the instrument to maintain the readings to within 


12 


TRANSMISSION 


-10 and 10 volts, the range of linearity. The resulting attenuation 
factor was 1.98 and its use or lack thereof is listed in the data file 
listing. 

3.23 The Data 

The data retrieval system used was a Campbell datalogger unit which 
stored the data in the format displayed in Table 3, where the 
parameters of each line reveal the date, the Z time in hours and 
minutes, seconds, and the output voltage reading, respectively. 



800 850 900 950 1000 

WAVENUMBER (cm-1) 


Figure 5: PRT-6 spectral response curve (Parsons 

deployment) 
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Table 3: 


Format of the PRT-6 data recorded in the 
datalogger. Although the format of the 
calibrated data is in a different form, both 
datasets can be accessed by the FORTRAN 
command "READ(*,*)". 


Day No. 

GMT 

Hrs., Mins. 

Sec. 

Voltage 

mv 

204 

2210 

52 

-1379 

204 

2210 

57 

-1383 

204 

2211 

2 

-1386 

204 

2211 

7 

-1389 

204 

2211 

12 

-1391 

204 

2211 

17 

-1392 

204 

2211 

22 

-1395 

204 

2211 

27 

-1397 

204 

2211 

32 

-1401 

204 

2211 

37 

-1404 


The calibration efforts of the PRT-6 output voltages into radiance 
measurements has been twofold. One method involved the 
intercomparison of the data with corresponding interferometer data 
taken at the same time and zenith angle and integrated over the 
aforementioned spectral transmission curve. The other method 
involved the acquisition of additional data taken from reference 
blackbodies. The combined results of both efforts have resulted in the 
curve of Figure 6. 








































PRT6 vs 897 - 935 Interferometer data 



Figure 6: Calibration curve. 


Raw voltage values may be reduced to radiance values by applying the following 
equation: 


I- (if) xo . 81418+646 . 80730 


where I is given in microwatts per square centimeter per steradian, A is the 
attenuation coefficient of 1.98 (used when applicable), and V is the output 
voltage signal in volts. With this in mind the following data files are 
available, shown here with their respective filenames, dates and times, zenith 
angle orientations, and electrical attenuation usage: 
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Filenames 


nol7p01 

nol7p02 

nol7p03 

nol8p01 

nol8p02 

nol8p03 

nol8p04 

nol8p05 

nol8p06 

no20p01 

no20p02 

no20p03 

no21p01 

no22p01 

no23p02 

no23p03 

no24p01 

no24p02 

no26p01 

no26p02 

no28p01 

no29p01 

de03p01, 

de03p02 

de04p01 

de04p02 

de04p03 

de05p01 

de05p02 

de05p03 

de06p01 

de06p02 

de06p03 

de07p01 


,dat nol7c01.dat 
,dat nol7c02.dat 
.dat nol7c03.dat 
dat nol8c01.dat 
dat nol8c02.dat 
.dat nol8c03.dat 
,dat nol8c04.dat 
.dat nol8c05.dat 
.dat nol8c06.dat 
.dat no20c01.dat 
.dat no20c02.dat 
.dat no20c03.dat 
.dat no21c01.dat 
.dat no22c01.dat 
.dat no23c02.dat 
.dat no23c03.dat 
.dat no24c01.dat 
.dat no24c02. dat 
.dat no26c01.dat 
.dat no26c02.dat 
.dat no28c01.dat 
.dat no29c01.dat 
dat de03c01.dat 
,dat de03c02.dat 
,dat de04c01.dat 
.dat de04c02.dat 
.dat de04c03.dat 
.dat de05c01.dat 
.dat d e0 5c02.dat 
.dat de05c03.dat 
.dat de06c01.dat 
.dat de06c02.dat 
.dat de06c03.dat 
.dat de07c01.dat 


[TCI 

End (UTC) 

Orientation 

m 

91 18:01 

17 Nov 91 19:27 

Vertical 

No 

91 19:32 

17 Nov 91 20:18 

Vertical 

No 

91 20:18 

17 Nov 91 21:19 

Vertical 

No 

91 19:23 

18 Nov 91 20:04 

Vertical 

No 

91 20:09 

18 Nov 91 20:27 

Vertical 

No 

91 20:37 

18 Nov 91 21:04 

Vertical 

No 

91 21:08 

18 Nov 91 21:33 

Vertical 

No 

91 22:00 

18 Nov 91 22:48 

Vertical 

No 

91 23:13 

18 Nov 91 23:44 

Vertical 

No 

91 22:10 

20 Nov 91 22:20 

Vertical 

No 

91 22:23 

20 Nov 91 22:47 

Variable 

No 

91 22:52 

20 Nov 91 23:55 

Variable 

No 

91 21:33 

22 Nov 91 04:06 

Vertical 

No 

91 14:42 

22 Nov 91 19:59 

Vertical 

No 

91 18:12 

23 Nov 91 18:43 

Vertical 

Yes 

91 19:00 

23 Nov 91 19:29 

Vertical 

Yes 

91 17:25 

24 Nov 91 22:44 

Variable 

Yes 

91 22:54 

25 Nov 91 02:43 

Vertical 

Yes 

91 14:20 

26 Nov 91 19:18 

Vertical 

Yes 

91 19:18 

26 Nov 91 23:45 

Vertical 

Yes 

91 15:42 

28 Nov 91 22:15 

Vertical 

Yes 

91 14:34 

29 Nov 91 17:33 

Vertical 

Yes 

91 19:50 

03 Dec 91 20:34 

Vertical 

Yes 

91 21:15 

03 Dec 91 22:05 

Variable 

Yes 

91 16:12 

04 Dec 91 20:29 

Vertical 

Yes 

91 20:29 

05 Dec 91 00:53 

Vertical 

Yes 

91 00:53 

05 Dec 91 05:07 

Vertical 

Yes 

91 00:53 

05 Dec 91 10:15 

Vertical 

Yes 

91 10:15 

05 Dec 91 15:35 

Vertical 

Yes 

91 15:35 

05 Dec 91 20:12 

Vertical 

Yes 

91 21:18 

06 Dec 91 03:51 

Vertical 

Yes 

91 07:39 

06 Dec 91 14:12 

Vertical 

Yes 

91 14:12 

06 Dec 91 20:25 

Vertical 

Yes 

91 15:41 

07 Dec 91 22:14 

Vertical 

Yes 

: for the uncalibrated and the 

calibrated data, 

and all 


dates and times are UTC. The variable zenith angle scans recorded on the third of 
December are unusual in that the respective sequence of zenith angles was 0, 15, 30, 5Q, 
60, and 75 degrees. Data files are archived at Colorado State University. 
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trueno.atmos.colostate.edu 

(129.82.107.109) 


and the following directory: 


/users/j ohnw/prt6 


3.3 Infrared Interferometer 



Figure 7: Infrared interferometer and calibration system. 

3.31 Instrument Description 

The instrument used in this experiment was a dual port emission 
interferometer manufactured by Bomem, Inc. It has an adjustable 
resolution ranging from 1 cm' 1 to 128 cm* 1 . FIRE data were taken 
with the interferometer at the 1 cm' 1 resolution setting. The detector 
is liquid nitrogen cooled. The useful range of the detector is from 500 
to 2000 cm' 1 . 

One port views a reference blackbody and is protected by a desiccant 
container. This reference blackbody is attached directly outside the 
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second interferometer port and it consists of a thermoelectrically- 
heated copper element which is grooved and coated to improve the 
emissivity. The reference blackbody temperature is monitored by a 
small thermocouple which is located within ten thousandths of an 
inch beneath the surface of the coated side. A simple toggling and 
monitoring control loop on the thermoelectric element maintain the 
temperature of the reference blackbody to within 0.1 C. 

Facing the target port of the interferometer is a rotatable gold-plated 
mirror mounted at a 45 degree angle with respect to the axis of the 
entering optical path. The position of this mirror is maintained by a 
computer controlled DC servo motor and control circuit, enabling the 
interferometer to scan the atmosphere at various zenith angles. 

The interferometer views the atmospheric emission through a KRS5 
(thallium bromide) window following reflection by rotatable gold- 
plated mirror. The beam then enters into a Michelson interferometer. 
The radiation from the reference blackbody in the reference port also 
enters the Michelson inte rferometer. This beam however is 180 
degrees out of phase with the atmospheric emission beam which 
results in the generation of a difference interferogram at the detector 
module. 
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3.32 Calibration and Measurement 

3.321 Blackbody sources 

Two blackbody sources are used to calibrate the 
interferometer. A cold blackbody source is a dewar of liquid 
nitrogen which is known to be at 77 K. The other blackbody 
used in the calibration is a warm blackbody source that is 
made out of copper. Its cavity is grooved and coated to insure 
a high surface emissivity. The copper blackbody is at the 
ambient air temperature; two small thermocouples located near 
the inner walls of the copper blackbody monitor the 
temperature. For the calibration of the instrument, it is 
necessary to record the emission at each wavenumber from 
each of the two blackbody sources. It is practical to use the 
liquid nitrogen since it is preferable that the two sources 
bracket the emission temperature of the atmospheric beam. 

The mirror is controlled by PC software such that the 
interferometer views the liquid nitrogen followed by the copper 
blackbody and then the atmosphere in sequence. 

3.322 Data Storage 

One data collection cycle includes viewing the two blackbody 
sources and the atmosphere. In the FIRE experiment each 
cycle consisted of 20 interferometric scans for each mirror 
position which are co-added to reduce the noise of the 
measurement. This yields three interferograms. The 
temperatures of the reference blackbody, the warm and cold 
blackbody sources, and the ambient air are measured every 0.2 
seconds. Every 100 temperature measurements are averaged 
and are then recorded at 20 second intervals on a CR21X 
Campbell datalogger. The temperature data are then 
downloaded to the PC which operates the interferometer. Each 
cycle takes slightly less than two minutes to complete, 
resulting in five temperature readings per object. 

3.323 Brief Outline of Calibration Procedure 

The data which are taken during each set of 20 scans by the 
interferometer is an interferogram. This may be written as: 
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Fix) - Q, exp[i<J>(v)] exp[i2nvx] cfv 


( 1 ) 


where: 


This interferogram can be processed by Fourier transform to a 
difference spectrum. The uncalibrated spectrum in terms of 
radiances may be written as: 


C, - |F| - r v (Iv + l*°) 


( 2 ) 


x is the optical path delay 

v is the wavenumber in cm' 1 

<p(v) is the phase response of the instrument 

r ¥ is the responsivity of the instrument 

L ¥ is the spectral emission 

L°„ is the offset emission of the instrument 

p is the complex Fourier transform 

Equation (2) shows the linear relationship between the 
uncalibrated spectrum and the spectral emission. The response 
of the detector and the offset are two unknowns to be derived 
from the calibration. These unknowns can be found from 
observing the cold and hot blackbody sources. 

The responsivity of the interferometer may be written as: 
r, - (3) 


and the offset required may be written as: 

14 - Sjsl£h-b V} W 

r. 


where B ¥ (T) represents the Planck function, the subscripts h 
and c refer to the hot and cold blackbodies respectively, and r 
refers to the reference blackbody. Now, with the known 
detector response and offset required, the atmospheric emission 
can be found and may be written as: 
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spectral radiance (W m-2 sr-l/cm-l) 


( 5 ) 



The interferograms are actually difference interferograms 
between the desired target and die reference blackbody, and 
since the reference blackbody temperature can be maintained 
to 0.1 C, the uncalibrated reference spectra (C^) subtracts out 
of equations (3) and (5). The calibration of the interferometer 
within the range of interest is within 2%. 

Interferogram and calibrated spectra data are archived at 
Colorado State University in a compressed format on exabyte 
tape. 


26 NOVEMBER 1991 at 1 600Z 



wavenumber [cm-l] 


Figure 9a: Example of a clear sky infrared spectrum observed on Nov. 26, 

1991. 
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pcetral radiance {W 


26 NOVEMBER 1991 at 2030Z 



wavenumber [cm - 1 ] 


Figure 9b: Example of an infrared spectrum observed in the presence of 

cirrus clouds on Nov. 26, 1991. 


References: 

Revercomb, H. E., et al., 1988: Radiometric calibration of IR Fourier transform 
spectrometers: solution to a problem with the high-resolution interferometer 
sounder. Applied Optics, 15, 3210-3218. 
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4.0 Solar Radiation 

(See Figure 3, section 3.1) 


4.1 Pyranometer 

4.11 Instrument Description 

The Eppley Precision Spectra Pyranometer (PSP) measures shortwave 
radiation including the direct component of sunlight, the diffuse 
component, and reflection from natural surfaces. 

It uses a multi junction thermopile blackened with Parsons black 
lacquer which is temperature compensated to give a response 
independent of ambient temperature. 

The pyranometers used had a pair of concentric hemispheres of Schott 
optical glass precision ground and polished. The inner hemisphere on 
both pyranometers was WG7 clear glass transparent from .28 to 2.8 
H m. On one pyranometer the outer dome was also WG7; the other 
used RG695 dark red glass transparent from .695 to 2.8 /im. 

The accuracy of the PSP is typically 1%. It should be noted, 
calibrations could have a diurnal and annual cycle and may change 
with temperature. 


4.12 Specifications 
Sensitivity 
Impedence 

Temperature Dependence - 

Linearity 
Response Time 
Cosine Response 


9 /xV/watt/m 2 
650 ohms 

± 1% over ambient 

temperature (20 * C to 
40*C) 

± .5% (0 to 2800 watt/m 2 ) 
1 second 

± 1% from normalization 
(0 to 70* zenith) 
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4.13 Data Reduction 


The instruments used in these experiments were compared with 
secondary standards to establish calibration coefficients relating the 
voltage output of the instruments to the irradiance. The following 
calibration coefficients were obtained and compared with pyranometer 
calibrations at NOAA CMDL: 


H(w/m 2 ) = aV + b 

Slope 

a(w m 2 mV 1 ) 

Y intercept 
b(w m‘ 2 ) 

Station 1 - WG7 

117.246 

-1.599 

Station 1 - RG695 

102.1832 

-1.726 

Station 2 - WG7 

100.9247 

-1.810 

Station 2 - RG695 

117.4345 

-1.659 


1 

4 


where V is the measured voltage output expressed in mv. 

For additional information concerning the data gathering system and 
data storage see section 3. 

Reference: 

The Eppley Laboratory: Instrumentation for the Measurement of the Components 
of Solar and Terrestrial Radiation. - 



Figure 10: Filtered and unfiltered pyrheliometer (see also figure 2, section 

2 . 2 ). 
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4.21 Instrument Description 


The Eppley pyrheliometer measures the direct component of solar 
radiation along a surface perpendicular to solar beam. 

The receiver is in the form of a thin 9 mm silver disc coated with 
Parsons optical black lacquer. Fifteen junctions of fine bismuth silver 
thermocouples are in thermal contact with, but electrically insulated 
from the lower surface of the disc. The cold junctions are in thermal 
contact with the copper tube of the instrument. A rotatable frame for 
three filters is provided. The receiver rotates allowing timed readings 
from each filter. The following filters were used: 


Channel 1 Schott filter 

Channel 2 Schott filter 

Channel 3 Schott filter 

Channel 4 Schott filter 

.53 - 2.8 fi m. 
.695 - 2.8 /ini. 
1.0 - 2.8 /xm. 

.3 - 2.8 fxm. 

Specifications 


Temperature Range 

-20 to 40* C (1% accuracy 
is achieved by a special 
temperature compensation 
circuit. 

Response Time 

98% of maximum output is 
achieved in 20 seconds 

Output 

3.0-3.5 /iv/cal/cm 2 /min. 

Internal Resistance 

400 ohms at 25 * C 

Emf 

8.8 2 x 10 -6 volts/watt/m 2 
(after comparisons with 
Eppley group reference 
standards under radiation 
intensities from 700 to 888 
watt/m 2 

Linearity 

within .5% of the above 
intensity 

Field of view 

5* 43’ 
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4.23 Data Reduction 


Two pyrheliometers tracked the sun using a LICOR 2020 solar tracker. 
Pyrheliometer data were reduced using a linear relationship between 
voltage output and radiance. This reduction equation is of the form 
N(w/m 2 ) = a * V (volts) where the constant a is 8.82E-06 w/m' 2 
Volt 1 and 6.905E-06 w/m' 2 Volt' 1 for pyrheliometers 1 and 2 
respectively. The three filters were installed in pyrheliometer 1 ; the 
second pyrheliometer made observations in the channel 4 spectral 
bandpass only. 

Data are archived at NCDS. 

References: 

The Eppley Laboratory: Instructions for the Installation of the Solar Tracker and 
Normal Incidence Pyrheliometer. 

The Eppley Laboratory: Standardization of Eppley Pyrheliometer (Normal 

Incidence). 

4.3 Multi Field of View Radiometer (MFOV) 



Figure 11: MFOV Radiometer 
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4.31 Instrument Description 

From the MFOV measurements, one may infer the optical depth of 
clouds and gain some insight into the equivalent MIE size 
characteristics of the scattering medium. This is achieved by 
measuring the solar irradiance in five different fields of view centered 
on the solar disk. 

The MFOV utilizes five silicon photodiodes (PV100A series) aligned 
along a common optical axis each collimated a different length to 
achieve five different full angle fields of view with a spectral bandpass 
.35 to 1.15 nm. 

Field of View Length fmm) Aperature Diameter (mm) 

(degrees) 


2' 

144.8 

2.5 

5’ 

57.2 

2.5 

10* 

29.2 

2.5 

20* 

14.0 

2.5 

28* 

10.2 

2.5 

4.32 

Photodiode Specifications 



The interior of the collimator is flat black to reduce scattering. The 
following are the characteristics of the PV100A photodiodes: 


2 

Active area - 5.1 mm 

Window thickness - 1.3 mm 

Window diameter - 6.1 mm 

Diode length (minus leads) - 4.2 mm 

Diode length (plus leads) - < 24.2 mm 


The MFOV tracks the sun with a LICOR microprocessor based solar 
tracker. 

The accuracy of the measurements are dependent on instrument 
calibration. Calibrations were performed at CSU on October 23, 1991 
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for clear sky cases. Data during the experiments are used to confirm 
this calibration data. 

For any additional information concerning the data gathering system, 
data storage, or data location. 

Data are archieved at CSU on exabyte tape. 


Reference: 

Raschke, R. A. and S. K. Cox, 1982: The determination of cloud optical depth from 
multiple field of view pyrheliometric measurements. Department of 
Atmospheric Science, CSU, Paper No. 361, 1-14. 

4.4 Sun Photometer 



Figure 12: Tracking Sun Photometer 

A MAINZ sunphotometer was deployed at the Parsons site in order to 
infer the optical thickness of cirrus clouds in the mid-visible portion 
of the solar spectrum. These data will also be used to provide 
calibration of the NASA lidar system. The photometer is not a 
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completely weather proof instrument and was deployed only during 
times of interest (see table 2) and only when precipitation posed no 
threat to instrument integrity. 

4.41 Instrument Description 

The MAINZ sunphotometer, model MZ II, is a precision electro- 
optical device which was originally designed to estimate 
aerosol optical depth. It was designed to be hand held or 
tripod mounted, and its output signal can be read from an LCD 
display on the top panel of the instrument. Several features of 
the instrument pertain to a "stand alone" type of measurement. 
In its deployment at Parsons the output of the instrument was 
obtained from a 0.0 - 2.0 volt output jack on the instrument, 
and several of the features (such as a selectable signal 
amplification, elevation angle measurements and peak signal 
holding feature) were not pertinent to its operation and will 
not be discussed here. 

The detector employed in the sunphotometer is an EG&G UV 
100 BQ, which is a UV enhanced photodiode linear to within 
1% for a range of incident irradiances spanning 7 decades. 
The dark current of the photodiode is 3 nA. 

A diopter with a piano convex lens and a surface mirror images 
the sun on an opaque screen in order to assure proper 
alignment of the instrument. 

A filter wheel allows the selection of 5 interference filters, 
although only the 0.500 ixm filter was used for this 
deployment. Table 4 summarizes the instrument capabilities. 
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Table 4. Technical specifications of the MAINZ MZ II - 85 sunphotometer 


Model 

MAINZ MZ II - 85 

Number of channels 

5 

Peak wavelength in /xm 

0.38, 0.412, 0.50, 0.675,0.862 

Wavelength accuracy 

±0.001 /xm 

Bandwidth 

0.002 - 0.005 pm 

Transmission 

15 - 50 % 

“ 

Full field of view 



Slope angle 

0.8- 

Filter diameter 

10 mm 

Effective filter diameter 

6 mm 

Type of detector 

EG&G, UV 100 BQ 

Relative signal accuracy 

0.5 % 


4.42 Instrument Deployment 


The sunphotometer was deployed on a LICOR model 2020 solar 
tracker. Actually, the instrument was originally deployed in tandem 
with the MFOV (Multiple Field Of View) photometer on 18 Nov 1991 
at 1900 UTC. Due to problems in keeping both instruments aligned, 
the instruments were placed on separate tracker units on 21 Nov at 
1400 UTC. The sun tracker was mounted approximately 1 meter 
above the ground and 10 m to the west of the H-Frame station along 
a berm which bordered the south shore of the power plant lake. 

As mentioned above, the sunphotometer is not designed for 
deployment in all types of weather and particularly not during 
precipitation events. When conditions warranted the instrument Was 
enclosed in a plastic bag containing a large supply of desiccant. The 
instrument was also covered in this manner during the night to 
protect it from moisture. Care was taken to keep the instrument 
aligned and viewing the solar disk; however, because of the 
photometers small field of view and imperfect alignment of the solar 
tracker the user of the data should beware of possible inaccuracies 
introduced by slight misalignment. Table 5 indicates the periods of 
data collection. 
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4.43 Data Collection and Archiving 

The data were collected by accessing a 0 - 2 Volt output jack on the 
side of the photometer. The data were recorded once every 2 
minutes, except as noted in Table 5. Data were recorded using a 
Campbell CR21X data logger. The data were downloaded to a DOS 
disk and converted to ASCII format. The data are archived in this 
original format and also in ASCII format on a UNIX workstation. 
Table 5 details the times of data collection as determined from an 
initial inspection of the archived data. All times are UTC. 


Table 5. Dates and times (UTC) of data collected by the Mainz Sun photometer. 


Date 

Start Time 

End Time 

Period 

21 Nov 91 

14:00 

15:35 

2.0 min 

25 Nov 91 

20:25 

22:00 

30 sec 

26 Nov 91 

14:00 

21:00 

30 sec 

03 Dec 91 

16:18 

16:48 

30 sec 

04 Dec 91 

15:16 

16:12 

30 sec 

04 Dec 91 

20:48 

22:48 

30 sec 

05 Dec 91 

20:24 

21:54 

30 sec 

06 Dec 91 

18:24 

21:18 

30 sec 


4.44 Calibration 


The data from the sunphotometer are in units of volts. No calibration 
exists at the current time for converting the signal to irradiance or to 
optical depth. Some information is contained in the data for clear sky 
conditions; however, further processing of these data is required to 
obtain calibration constants. A laboratory calibration will be 
performed on the sunphotometer using an irradiance standard lamp 
during spring of 1992. 
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5.0 Upper Air 


5.1 Radiotheodolite 

5.11 Instrument Description 

At the Parsons site, rawinsondes were launched primarily in support 
of surface radiation measurements. The sondes were Intellesondes, 
made by A.I.R. (Atmospheric Instrumentation Research, Inc.). The 
sondes transmit at 1680 MHz, and contain a capacitance aeroid 
barometer, a 2mm diameter rod thermistor, and a carbon hygristor. 
The barometer has a range from 1050 to 5 mb with a sensor precision 
of 1 mb. The range of the thermistor is from 50 * C to -90 * C with 
a precision of 0.5* C. The carbon hygristor has a range from 5% to 
100% RH with a precision of 3% RH. (Note that a humidity of 20% 
is often observed as a minimum value, especially at low pressures and 
cold temperatures.) 

The sonde radio signal was tracked by an A.I.R. automatic 
radiotheodolite (model AIR-3 A-RT). The antenna consists of an eight 
element phased array. From the antenna azimuth and elevation 
angles the winds can be calculated. The antenna has a static tracking 
accuracy of 0.5 degrees. (Due to hardware problems, the tracking 
accuracy varied and was generally less reliable than the 
specifications.) 

5.12 Data Format/Ascent Schedule 

The rawinsonde data consist of julian day, hour (UTC), minute, 
second, height (km), pressure (mb), temperature (K), relative 
humidity (%), azimuth angle, and elevation angle, and is stored with 
the format (3I3,F6.2,6F10.3). An example of the data is seen in 
Table 6. Table 7 shows the sonde flight times and brief comments on 
the data quality. The data are available from NCDS (NASA Climate 
Data System). 
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Table 6. 


Sample of a rawinsonde data file. 


JD HR MIN SEC HT PRES(mb) TEMP('A) Relative Azimuth Elev. 

(KM) Humidity % deg. deg. 


330 

14 

58 

20.08 

0.269 

988.320 

278.620 

67.090 

337.830 

- 0.730 

330 

14 

58 

22.23 

0.269 

988.330 

278.620 

68.140 

338.580 

- 0.730 

330 

14 

58 

26.46 

0.269 

988.300 

278.600 

67.340 

337.920 

- 0.520 

330 

14 

58 

28.60 

0.269 

988.340 

278.580 

67.070 

337.730 

- 0.930 

330 

14 

58 

30.68 

0.269 

988.330 

278.580 

67.410 

338.000 

- 0.680 

330 

14 

58 

32.83 

0.269 

988.330 

278.580 

67.090 

338.270 

- 0.680 

330 

14 

58 

34.91 

0.269 

988.270 

278.610 

67.210 

338.350 

- 0.700 

330 

14 

58 

37.06 

0.269 

988.310 

278.630 

66.940 

338.520 

- 0.550 

330 

14 

58 

39.14 

0.269 

988.340 

278.640 

66.690 

338.550 

- 0.730 

330 

14 

58 

51.83 

0.340 

979.720 

278.050 

64.260 

341.850 

- 1.020 

330 

14 

58 

53.97 

0.356 

977.870 

277.910 

64.460 

334.770 

2.400 

330 

14 

58 

56.06 

0.371 

976.060 

277.770 

64.690 

325.380 

10.520 

330 

14 

59 

0.29 

0.402 

972.350 

277.620 

65.200 

324.700 

15.900 

330 

14 

59 

2.38 

0.416 

970.630 

277.560 

65.310 

325.000 

21.300 

330 

14 

59 

4.52 

0.432 

968.680 

277.430 

65.450 

330.500 

21.330 

330 

14 

59 

6.61 

0.446 

966.990 

277.310 

65.670 

339.520 

21.330 

330 

14 

59 

8.75 

0.465 

964.720 

277.260 

66.680 

341.300 

21.330 

330 

14 

59 

17.21 

0.524 

957.820 

277.320 

68.260 

348.950 

31.150 

330 

14 

59 

19.35 

0.536 

956.390 

277.230 

68.430 

349.380 

30.420 

330 

14 

59 

21.49 

0.552 

954.520 

277.200 

69.090 

350.580 

29.850 
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Table 7. Rawinsonde Launches 


Date Time Lowest Comments 

(UTC) Pressure 


911113 

911113 

911117 

911118 

911118 

911118 

911118 

911119 

911119 

911120 

911120 

911120 

911121 

911121 

911121 

911121 

911122 

911122 

911122 

911122 

911124 

911124 

911125 

911125 

911125 

911125 

911126 

911126 

911126 

911128 

911129 

911203 

911204 

911204 

911204 

911205 


18:28 

21:20 

18:48 

15:45 

18:58 

21:40 

22:41 

15:09 

17:33 

19:03 

21:46 

23:01 

16:38 

18:16 

21:06 

22:52 

02:00 

13:40 

15:48 

18:11 

18:15 

21:55 

00:09 

00:45 

15:12 

20:08 

14:58 

17:43 

19:36 

19:29 

16:22 

21:06 

16:54 

20:41 

22:15 

03:55 


250 mb 
207 
26 
25 
28 

125 

205 

963 

193 

83 

131 

46 

407 

30 

190 

222 

61 

109 

80 

276 

39 

44 

734 

46 

66 

78 

124 

94 

63 

82 

39 

341 

516 

597 

155 

384 


No wind, data missing between 871 and 783 mb 
No wind 


Data missing 915 - 184 mb (data collection 
problem) 

Data missing between 407 and 163 mb 

Data collection problem 

No wind, data missing 325 -273, 263 - 193 mb 

No wind 

No wind 

No wind 

No wind, data missing 621 - 587, 587 - 429 mb 


No wind 
No wind 
No wind 
No wind 
No wind 
No wind 
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Table 7, Continued 


Date 

Time 

(UTC) 

Lowest 

Pressure 

Comments 

911205 

14:25 

197 


911205 

16:17 

49 


911205 

18:23 

29 


911205 

21:07 

141 


911206 

01:39 

46 


911206 

05:29 

615 


911206 

06:45 

52 


911206 

12:30 

38 

No wind, data missing between 212 and 174 mb 

911206 

17:04 

25 

No wind 


Note: Many times immediately after launch the theodolite initially had trouble tracking 
causing some low level erroneous winds. 


5.2 Model 400 Wind Profiler 



Typical beam configuration in wind profiling consists 
of three beams: one vertical, and two tilted IS * 
from the zenith (to the east and north, for example). 
Under some circumstances, two additional beams 
are needed (such as to the south and west ). 


Figure 13: Wind Profiler System 


CHAIN UNK PENCE 85 X 65 120 X JO m*wii 


ANTENNA GROUND PLANE *6' X 46 (14 X 14 rrwtlKI} 


ANTENNA ARRAY 39' da (12 malar da ) 



I' 
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5.21 System Introduction 

The model 400 wind profiler provides measurements of wind 
(horizontal and vertical components) above the radar site up to 1 5km. 
The whole system mainly consists of two parts: the radar system and 
the computer work station for system control and data processing. 
The data are provided continuously every 6 or 10 minutes (set by 
operators). Data can be used on-site, saved on magnetic media and 
transferred to some central system. The accessed data processing 
system provides the following format data: 

5.211 Spectrum every 6 or 10 minutes (Figure 14). 

5.212 Processed data including signal power, velocity, and signal 
spectral bandwidth every 6 or 10 minutes (Figure 15). 

5.213 Hourly and 6-minute or 10-minute horizontal wind profile 
(Figures 15 and 16). 

5.22 System Principal 

The wind profiler is a Doppler radar which depends upon the 
scattering of electromagnetic energy by small irreg ularit ies in the 
refractive index. The spatial variation in this refractive index 
encountered by propagating microwave energy causes a small amount 
of the energy to be scattered (or dispersed) in all direction of which 
the backscattering is received by the radar. The irregularities, which 
have the size comparable to one-half the radar wavelength, are the 
most effective sources of scattering. Most of the energy propagates 
through the refractive irregularities without being scattered. 

The transmitted microwave energy is at a frequency of 404.37MHz. 
The motion of the air toward or away from radar will cause a 
frequency shift of the received signal. T his frequency shift is 
proportional to the speed of the air motion, and can be used to 
calculate the wind speed. Following is the Tycho algorithm used to 
calculate horizontal wind speed: 

V re = U sin 75 0 + w cos 75* 

V m = V sin 75' + w cos 75* 

V ra = w 

Where, V re , V m> and V n are radial velocities measured in the east, 
north, * and zenith directions, respectively. U, V are east-west, 
north-south components of horizontal wind, respectively; W is the 
vertical velocity. 
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Typical power spectrum showing measurement parameters. 


Figure 14: Spectrum 



VUDCm (p/a) 


An example of wind data for one six-minute integration. 
Profiles of the three components (u, v. w) are shown separately. 
Also shown is the horizontal wind speed , which is calculated 
from the east-west (u) and north-south (v) components. 


Figure 15: Ten minute wind profile 
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An example of hourly averages of horizontal winds . 
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Figure 16: Hourly averaged wind profile 

5.23 Specifications 

5.231 Wind Profiler Site Location: 

Parson, Kansas: Elevation: 269m Latitude:37* 18’ 

Longitude:95 * 0 T 

5.232 The Wind Profiler Parameters: 

Table 8 lists the operator selected parameters which determine 
the operating characteristics of the wind profiler. Values 
assigned to these parameters are given in the Appendix. Here, 
we should mention that, during the whole FIRE period, we 
used five beams. So the data time cycle was about 10 minutes 
plus one minute to transfer the spectrum to an optical disk 
archive. 

5.233 Height Coverage and Resolution: 

The wind profiler has high and low modes. The low mode 
covers the height from 500m above ground level (AGL) to 
9.25km AGL with height resolution of 250m. The high mode 
covers the height from 7.50km AGL to about 1 6.25km AGL 
with a height resolution of 1000m. 
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Table 8. 


List of the operator selected parameters which determine the operating 
characteristics of the wind profiler. 


PARAMETERS 

SHORT DESCRIPTION 

ANTDIR 

Antenna direction ? 

ARA 

Anti-range aliasing for each mode 

CCP_GAIN 

Digital gain of the real and quadrature channels with the A/D 
converter 

CODE 

Complementary coding for each mode 

CON_DELTAV 

Consensus velocity half-width (m/s) 

DELAY 

Delay from fall of TX pulse to first range gate (1/6 /is) 

DODC 

Enables DC removal 

DOPOWCOR 

Corrects signal power values due to the range dependent 
attenuation of the T/R switch 

FILTER 

Receiver filter for each mode 

GCNUM 

Number of points around zero used in ground clutter removal 
(same for all modes) 

POW_FACTOR 

Power scaling factor 

PRP 

Pulse Repetition Period (0.1 /is) 

PW_BAUD 

Pulse width for each mode (1/6 /is) 

RXDELAY 

Receiver delay (0.01 /is) l 

RXGAIN 

Receiver gain 

SAMPLE_MIN 

Minimum number of agreeing samples for consensus 

SPACE 

Time between range gates for each mode (1/6 /is) 

SPAN 

Number of spectral points around zero that are searched for 
peak power 

SPECTRA 

Number of spectra averaged for each data-taking mode 

SSTDODC 

Enables SST DC-removal 

TDA 

Number of samples per time domain average for each mode 

WINDOW 

Enables a HANN windowing function for the FFT 
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5.234 Horizontal Spatial Resolution: 

To use Tycho wind calculation algorithms, it is assumed the 
flow to be horizontally homogeneous between beams. 
Therefore one can only calculate the wind field which has 
spatial scale at least bigger than the distance between north 
(south) and east(west) beam volumes. This distance can be 
refereed as the horizontal spatial resolution when we neglect 
the air motion. At 10km AGL, this distance is about 3.8km. 

5.235 Direction Deviation: 

When we use above Tycho algorithms, we assume the wind 
profiler is aligned along the N-S direction, actually, there was 
a 4.5° deviation east to north. Therefore the true wind 
component U’ and V are given by: 

U’ = U sin 4.5’ - V cos 4 .5 * 

V = U cos 4.5* + V sin 4.5’ 

5.236 Wind Profiler Operation Times: 

The system operated from 12:00 UTC on Nov. 13, 1991 to 
Dec. 7, 1991; at times within this period, mode 6 was used to 
measure temperature by RASS once each hour so that the wind 
profile data of this mode was lost. Also, there were 
scheduled-off times when satellites passed over the site. 

5.24 Data Description 

5.241 Archiving of the Data: 

The data were archived and the associated format and media 
used are given below. 

5.2411 Raw spectral da ta: 

Raw spectral data were archived on optical disks 
(following time domain average, FFT and 
spectral domain average, but before other signal 
processing) for each gate and each cycle. The 
file names are in the format: yymmdd_hhmmss. 
***. Here, yy, mm, dd, hh, mm and ss mean 
year, month, day, hour, minute and second, 
respectively. The data format in each file is: 
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Title: 

mode 1 _26_N o v_l 991_1 2:32 :44_Gate_24_Height_xxxxx/ 

Byte 1 --4: ’mode’ 

Byte6--7:mode number 1--10 

Byte9,10:day 

Bytel2-14:month 

Bytel6-19:year 

Byte2 1-28: time 

Byte30-33:’Gate’ 

Byte35-36:gate number 1--36 
Byte38-43:’Height’ 

Byte45 -49: height in digital format 
Byte50:retum 

After above title, a 256x4 byte binary array is saved. This array is the 
spectrum for the mode and gate specified by the title. In each file, we have 
360 arrays. 

5.2412 Processed data 

The processed data are derived from the Tycho processing 
program. It includes all of the moments data: signed power, 
radial velocity and speed variance. The data have been saved 
in the TK50 tapes with the following file format: 

yymmdd_hhmmss.new (covers the every cycle data in one hour) 
yymmdd_hhmmss.pro (covers the data which has been hourly 
averaged) 

and, 

yymmddhh.ard (cycle data after quality control) 
yymmddhh.arc (hourly-averaged data after quality control) 

5.242 Accessing of the Data 

As mentioned above, spectra data are archived at CSU on 
optical disk and moment data on TK 50 and exabyte tapes. 
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van de Kamp, D. W., 1988: Wind Profiler Training Manual #1 -Principles of Wind 
Profiler Operation. NOAA Documents Developed for the National Weather 
Service Office of Meteorology. 
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5.3 Radio Acoustic Sounding System 
5.31 Instrument Description 

The Radio Acoustic Sounding System (RASS) is used to acquire 
temperature profiles of the boundary layer. There are several 
methods available to use RASS all of which involve introducing an 
audio signal into the atmosphere and tracking the speed of the 
acoustic wavefronts using a sensitive Doppler RADAR. By adjusting 
the audio frequency such that the acoustic wavelength is half that of 
the transmitted radio frequency (i.e. satisfies the Bragg condition) and 
by using a wind profile RADAR (or wind profiler), the temperature at 
a given height can be obtained from the equation (Nalbandyan, O. G., 
1977) 


C a -20.047/7>v w (mis) 


where C a is the speed of sound, T v is virtual temperature in degrees 
Kelvin and v w is the wind velocity component parallel to acoustic 
wave propagation. 

The RASS method used presently at CSU is the Receiver Offset 
Method. After a user-determined number of profiler cycles, the 
intermediate frequency in the wind profiler's receiver is offset via a 
single-sideband up-converter, (Strauch, R. G., 1988). The offset 
causes the baseband of the receiver's discrete-time filter to be shifted 
to the audio frequency used by RASS. 

Shortly after the offset is applied and before the RASS mode begins, 
a separate oscillator sweeps through a narrow band of audio 
frequencies selected according to expected air temperatures. This 
signal is then amplified and broadcast over four loudspeakers which 
are placed approximately 10 meters away from and encircling the 
RADAR antenna. When the Bragg condition is satisfied, the strongest 
signals at the receiver for each range gate will be those of acoustic 
wavefronts traveling at the velocities corresponding _ to the 
temperatures at those altitudes. The Receiver Offset Method was the 
RASS procedure employed during the FIRE II experiment in Parsons, 
KS. 
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5.32 Specifications 

5.321 Electromagnetic 


Radar 

Tycho DORA 400S Wind 
Profiler 

Transmit Frequency 

404.37 MHz 

Radiated Power 

40 kW (peak) 

Beamwidth 

5 degrees 

Antenna Type 

Phased Array Coaxial 
Collinear Dipoles 

Antenna Diameter 

12 meters 


5.322 Acoustic 

Transducer Type 

Transducer Sensitivity 
Radiated Power 
Beamwidth 


Folded Horn with 
Compression Driver 

105 dB SPL (lW/lm) 

75 Watts (average) 

110 degrees horizontal 

130 degrees vertical 


5.323 Data Summary 


All data collected by the wind profiler are stored in two 
formats - as processed spectral moment data and as raw 
spectra. Temperature profiles are determined from the spectra 
(see Figure 17) by the equation 



20.047 ) 


(ft 


where T v is virtual temperature and v is the apparent 
velocity reported by the profiler. The offset velocity, v ofh> 
results from the offset frequency, f offs , and is calculated from 
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v +--yf<* Ms) 


where A em is the wavelength of the electromagnetic radiation. 


During the FIRE II experiment, RASS data were collected every 
sixth profiler cycle during normal daylight operating periods. 
During intensive data collection intervals or during 
maintenance, RASS may have been activated more or less 
often. 
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Figure 17: Sample RASS spectra (colder temperatures to the right). 
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6.0 Goud Measurements 


6.1 Laser Ceilometer 

6.11 Instrument Description 

The ceilometer utilized at the Parsons site of the 1991 FIRE IFO 
experiment is an active LIDAR system which employs a vertically 
oriented Gallium-Arsenide laser to detect cloud base height. Relating 
specifically to the FIRE experiment, it was utilized to provide a record 
of when low-level cloudiness was present over the Parsons area. 

Within the time frame of the experiment it was not unusual for the 
atmosphere to present a situation in which a desired cirrus layer study 
was largely obscured by low-level cloudiness. It is in such situations 
that, often to the frustration of the observers, radiometric cirrus 
studies become obscured by such clouds even in situations of partial 
coverage. The ceilometer was deployed to give a continuous database 
describing when such low-level cloudiness was present. It 4 is also 
hypothesized that the ceilometer may have also occasionally detected 
an actual cirrus cloud near the top of the instrument’s vertical range. 


6.12 Specifications 

The ceilometer utilizes a 20 watt near-infr ared Gallium -Arsenide laser 
operating at a wavelength of .91 microns. It employs the use of 1024 
range gates which give it a vertical resolution of 25 feet up to a 
maximum range of 25,600 feet, the approximate lower edge of the 
vertical range of cirrus cloud ocurrences. The fields of view of the 
transmitter and receiver have been estimated to be of the order of 1 
degree. 

The temporal resolution can be set at one of several settings between 
30 seconds and 12 minutes, as the instrument runs through a cycle of 
of data collection, interpretation, and transmission. This begins with 
a collection mode where the ceilometer, taking a reading from all of 
the range gates 1000 times per second, first performs a noise level 
analysis by taking three such readings in which the laser does not fire. 
All these measurements are then condensed into one noise level 
reading. The ceilometer then takes 5120 of the same sets of readings 
from each range gate, this time with the laser pulsing. During this 
operation a histogram of the received signal is created in the 
following manner: 

The signal from each range gate is compared with the previous noise 
level reading, and depending on whether the received signal is above, 
below, or at the given noise level a value of 1, -1, or zero is, 


46 



respectively, added to the output histogram value for that specific 
range gate. This is performed for each of the 5120 pulsed readings, 
giving a final output in the form of a simple array of integers. 

After performing a peak-location algorithm which determines the 
height of the cloud base (if one exists), the data are then transmitted 
to a 286 desk top computer where it is stored and displayed in real 
time. The ceilometer then goes into a wait mode until sufficient time 
has elapsed for another reading to be done according to the set value 
for the temporal resolution. 

6.13 The Data 

The data can be stored in either binary or ASCII, and at present all of 
the data are available in either form. In its ASCII form the data takes 
the following format: 

-7 

-7 

-7 

0 

14 

36 

0 

6 

7 

12 

1991 

86 

-400 

-1 

-1 

-1 

178 

217 

211 

13 

Each file contains a series of output histograms, each containing 1050 
integers. The first 1018 entries are the histogram values for each 
height gate (25 feet) (meaning that only information up to 25,450 
feet is actually recorded), and the final 32 integers display the 
following information: 
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Position 

Function 

Example 

1019-1023 

* 


1024 

(Always zero 

0 

1025 

Hour 

7 

1026 

Minute 

32 

1027 

Second 

0 

1028 

Day of Week 

4 

1029 

Date 

30 

1030 

Month 

5 

1031 

Year 

1991 

1032 

Waveform Mode 

86 

1033 

Cloudbase 

-1 

1034-1036 

* 


1037 

Laser Temperature 

180 

1038 

Receiver Temperature 

260 

1039 

Ambient Temperature 

250 

1040-1045 

* 


1046 

RMS Noise 

408 

1047-1050 

* 



where a cloudbase value of -1 indicates clear sky, the temperatures are given 
in tenths of a degree Celsius, and an asterisk indicates dummy parameters. 

Due to inclement ambient temperatures within the ceilometer circuitry 
caused by a malfunctioning heater, there were apparent electronic problems 
within the unit resulting in abrupt "spikes" in the histogram data. A simple 
program was created which filtered out these unrealistic values, replacing 
them with data linearly interpolated from adjacent values. Also, since these 
data "spikes" have before been previously interpreted by the ceilometer as a 
false cloud base, whenever such filtering has occurred within a given 
histogram, the cloudbase parameter (1033) was changed to a value of -400 
(-10,000 ft). A record was also created of all such changes made for each 
file. ’ 
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The following data files are available: 


Filename 


Start fUTCi 


End fUTCI 


Resolutioi 


10nov91.dat 10 Nov 91 20:48 

llnov91.dat 11 Nov 91 15:36 

12nov91.dat 12 Nov 91 14:12 

13nov91.dat 13 Nov 91 14:12 

14nov91.dat 14 Nov 91 15:00 

15nov91.dat 15 Nov 91 14:48 

16nov91.dat 16 Nov 91 16:12 

17nov91.dat 17 Nov 91 16:12 

18nov91.dat 18 Nov 91 14:12 

19nov91.dat 19 Nov 91 14:48 

20nov91.dat 20 Nov 91 16:12 

21nov91.dat 21 Nov 91 15:00 

22nov91.dat 22 Nov 91 13:12 

23nov91.dat 23 Nov 91 17:36 

24nov91.dat 24 Nov 91 17:12 

24nov291.dat 25 Nov 91 01:12 

25novl91.dat 25 Nov 91 16:39 

25nov291.dat 25 Nov 91 23:09 

25nov391.dat 26 Nov 91 00:08 

26novl91.dat 26 Nov 91 15:16 

26nov291.dat 26 Nov 91 17:21 

26nov391.dat 26 Nov 91 21:36 

26nov491.dat 26 Nov 91 22:36 

27novl91.dat 27 Nov 91 19:00 

28novl91.dat 28 Nov 91 18:00 

28nov291.dat 28 Nov 91 19:10 

28nov391.dat 28 Nov 91 22:12 

29novl91.dat 29 Nov 91 14:48 

30novl91.dat 30 Nov 91 14:48 

01decl91.dat 01 Dec 91 17:00 

02decl91.dat 02 Dec 91 16:24 

03decl91.dat 03 Dec 91 17:12 

04decl91.dat 04 Dec 91 16:24 

05decl91.dat 05 Dec 91 15:12 

06decl91.dat 06 Dec 91 14:48 

07decl91.dat 07 Dec 91 14:36 


11 Nov 91 13:48 

12m 

12 Nov 91 14:00 

12m 

13 Nov 91 13:48 

12m 

14 Nov 91 14:48 

12m 

15 Nov 91 14:36 

12m 

16 Nov 91 16:00 

12m 

17 Nov 91 16:00 

12m 

18 Nov 91 14:00 

12m 

19 Nov 91 13:48 

12m 

20 Nov 91 15:48 

12m 

21 Nov 91 13:48 

12m 

22 Nov 91 13:00 

12m 

23 Nov 91 17:00 

12m 

24 Nov 91 15:24 

12m 

24 Nov 91 17:24 

12m 

25 Nov 91 14:40 

lm 

25 Nov 91 23:07 

lm 

26 Nov 91 00:02 

lm 

26 Nov 91 04:08 

4m 

26 Nov 91 17:08 

4m 

26 Nov 91 21:13 

lm 

26 Nov 91 22:31 

lm 

27 Nov 91 15:48 

4m 

28 Nov 91 15:00 

12m 

28 Nov 91 19:00 

12m 

28 Nov 91 22:08 

lm 

29 Nov 91 14:24 

12m 

30 Nov 91 14:36 

12m 

01 Dec 91 16:48 

12m 

02 Dec 91 16:12 

12m 

03 Dec 91 17:00 

12m 

04 Dec 91 16:12 

12m 

05 Dec 91 14:48 

12m 

06 Dec 91 14:36 

12m 

07 Dec 91 14:12 

12m 

08 Dec 91 16:12 

12m 


The files are available in the following directory: 

/users/johnw/ceil 

and appear with a first letter which is either "u", "f or "s", revealing files 
which respectively contain either unfiltered data, filtered data, or a summary 
of the filtering process. 
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6.2 Video All Sky Camera 



Figure 18: Video all sky camera convex mirror. 

In order to archive a visual record of sky conditions, CSU’s wind profiler site 
included a VCR and camera facing an upward looking convex parabolic 
mirror. The camera was mounted over the upward-facing mirror, thus giving 
a full-sky view— obscured only by the reflection of the camera. 

In all, 25 six-hour VHS video tapes were recorded over the course of the 
experiment. The camera was only operated during daylight hours when sky 
conditions were acceptable, i.e. not pouring rain. Each tape covers a six 
hour period. The video tapes are archived at CSU. 

6.3 NASA Langley 8” Cloud Lidar System 

(Operated by Joe Alvarez, NASA Langley) 

6.31 Instrument Description 

The lidar deployed by Langley Research Center (LaRC) at a field site 
near Parsons, Kansas for the FIRE IFO during November 13 to 
December 6, 1991 was the LaRC 8" Lidar System housed in a 
recreational vehicle. This lidar system is comprised of a frequency 
doubled Nd:YAG laser serving as the lidar transmitter and an 8" 
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Cassagrain depolarization sensitive receiver system attached to a 386 
based data acquisition system. 

6.32 Operating Parameters 

The table below presents the operating parameters of this lidar system 
as it was configured and used for the FIRE Kansas IFO. 


TRANSMITTER 


Wavelength 532 nm 

Pulse Energy 150 mJ 

Pulse Rate 10 Hz 

Pulse Width 6 ns 

Beam Divergence 0.6 mrad 

Beam Diameter 6 mm 


RECEIVER 


Telescope Diameter 
Field of View 
Filter Bandwidth 
Detector, PMT 
Quantum Efficiency 
Calibration 


20 cm 
2 mrad 

1 nm, 50% transmittance centered at 532 nm 

EMI 9658 

15% 

Zero order half-wave plate mounted on rotatable 
assembly used for depolarization calibration 


DATA ACQUISITION AND ANALYSIS SYSTEM 


Computer 
A/D Resolution 
Altitude Points 
Altitude Limits 
Altitude Resolution 
Storage Media 
Data Analysis 
Graphics Plotter 
Lidar Record 


386 system, 20 Mz 

10 bits 

4096 

0 to 60 km, 2 to 15 km normal data range 
30 m(amplifier bandwidth usually 5 MHz) 

Optical Disk, 400 Mb per disk 

486 system, 33 Mz 

286 system, 12 Mz and inkjet printer 

Elementary unit of recorded data consists of hardware 

average of 150 lidar pulse responses (15 second average) 


Data acquisition commenced at 6:00 am on Wednesday, November 13, 1991 and 
continued intermittently until 11:29 am Friday, December 6, 1991. Approximately 
.700 Mb of raw data were collected at the field site during a cumulative 143 hours 
of operation. As mentioned above, the basic unit of data collection were 150 lidar 
pulse responses averaged over 15 second intervals. 
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Lanqley Research Center 8" Lidar 

1 991 FIRE IFO, Parsons, Kansas 



Figure 19: Presents operating times of the LaRC 8" Lidar System while deployed at Parsons. The "other obs" refers 

to lidar operation occurring during times where no cirrus were detected. During these times either the 
sky was clear or the detected low level clouds were obscuring overlying cirrus. 
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which will shortly become available at the NCDSS. 
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RX-GAIN 
SAMPLE-NUN 


SPACE 
SPAN 
SPECTRA 
SSTDODC 
1 TDA 


WINDOW 
DO POW COR 


MODE 1 MODE C 


NEW 


MODES I MODE 4 II MODE 5 I MODE 6 


11223 


NEW I NORM I HEW ||rOUA| NEW 1 NOW | NEW 


|.MEgal 




KV 




>AW»i 


COMMOM WWM5 |UNtt 


BLANK -TX 
CONSESUS-SET 


COOLDOWN 

FREEDELAY 

NUMGATE5 

NUMMDDES 

•TR.TX 


TX-BIANK. 

TX-DUTY 

TX-TR 


B 


FMN 




























parameter 


ANT_ Did 
ARA 

CLP- GAIN 
CODE 

CON-DELTAV 


GaCWLiM 
ROW. FACTOR. 


PW-BAUD 
RX_ DELAY 
RX-GAIN 
SAMPLE-MIN 
SPACE. 

SPAN 

SPECTRA 

SSTDODC 


WINDOW 

DOPOWCjOR. 


N\OOE7 IVAODE B I NlDOLQ 


NKMSMR NEW |NOft£l 




MODE 10 I MODE 11 I AAODE IE 


ndehii mem || nob* I now h ndeh }wtw 


I 


im 

iiunm 























M 


MODE < 


uwrrl morim new 


PARAMETER 


ANT-DIE 
ARA 

CCP-GAIN 
CODE 

CDN_ DELTA V 
DELAY 300 

DO DC o/l 1 

FILTER cooe / 

GCNUM N 3 

RQW_ FACTOR. 45U Kfled 
"PRP 

PW_ BAUD 
RX- DELAY 
RX-GAIN 
SAMPLE-NUN 


SPACE 
SPAN 
SPECTRA 
SSTDODC 
TDA 


WINDOW 

DOPOWCOR 




DATE OF CHANGE: j^jgi 


MD0E3 1 MODE 4 


KAODF 5 1 W3DE6 


NEW g I M6W 1 


* < 

1 

0 • 


1 ^ 
1 


10 

63 <*v 


MJEN NEW 


m? 

208 

2Srf 


O 7 1 £3 




1331 


V'H - Is- Y ?r 






































PfcZAMETEK 


AUT-DIE 

ARA. 

CC.P-6AIN 

CODE 

CON. Daw 


cornu 


W0DE1 


uwrr I MOAftM NEW I NORM. 


DATE OF rHANfiF: m?o&.f*r 


NDEM. I NEW 


MODE 4 


MODE 5 I MODE-6 


rCKNVl NEW I NDCN I NEW 


2 



PRP 

PW.BAUD 
RX- DELAY 
RX-GA1N 
SWSAPLE-hMM 


SPACE 

SPAN 

SPECTRA 

SSTDODC 

TDA 


WINDOW 

DOPOWCOR 





lira 


BLANK _TX 
C0N5ESU1.SET 


COOLDOWN 
FREEDELAY 
NUNVG»ATES 
NUNMDDES 
TR.TX 


TX-BLANK 

TX-DLTTY 

TX-TR 


^ •- jsi ?r 



























parameter R UNIT 


ANT_ DIR. 

ARA 

CCP-GAIN 

CODE 

CO N-DELTAV 
DELAY 
DO DC 
FILTER 
GCNUNA 
ROW. FACTOR. 


PRP 

PW-BAUD 
I RX- DELAY 
RX-GA1N 
SAMPLE- MIN 
SPACE 
SPAN 
SPECTRA 
SSTDODC 
TDA 


WINDOW 

DOPOWCjOR. 


t^DDElI NAPPE 8 I NAPPE 9 1 MODE 10 1 NODE H IMOOEIE 


ndcmh NEW | norm 




ISOeKH NEW B NOB* | NON || NDR*|VjeW 


Jf-Wi I 





















MODE 7 I MODES I NODE 9 


(WOOL 10 1 MODS 11 | MOOt 12. 


NOCMR N£W II NORM 


ANT_ DIR. 

ARA 

CCP-GA1N 

CODE 

CON-DtLTAV 
DELAY 
DO DC 
FILTER 
GCMM 
POW- FACTOR 


PRP 

PW-BAUD 
RX_ DELAY 
RX-GAIN 
SAMPLE- MIN 
SPACE 
SPAN 
SPECTRA 
SSTDODC 

TDA 

’WINDOW 
DOPOWCjOR. 


NDCMlI MEW H NORM I NOW || NO KM |VJ€Mi 




| 


tea 






















EIECS 


rnu-mrT: £&SS«u.D«rE OF CHANGE: »A Ml 


MODE! 



MODE 4 II MODE 5 I WODf.6 


PARAMETER 


AKII_ DIR 
ARA 

CCP-GAIN 

CODE 

con.dqtav 


DELAY 
DO DC 
FILTER 
GCNUM 
FOW_FACTOR 


PRP 

PW.BAUD 
RX_ DELAY 
RX-GA1N 
SAMPLE- IAIN 


SPACE 

SPAN 

SPECTRA 

SSTDODC 

TDA 


WINDOW 
DO POW COR 




r«XNV I NEW I NDV* | NEW 


COMMON PAKM& |UNH 


BIANIC-TX 

fnMv^K-ser 


COOLDOWN 

FREEDELAY 

NUNMjATES 

NUNlfWDDtS 

TR.TX 


TX- BLANK 

TX-DUTY 

TX_TR 




10 \ 

(3V & 
281 3>S 

*1 60 































parameter 


i^ocei tvAQpe a n\qdc Q I mode. <o 1 mope n moos <s 

UNIT II NDKMl NCI*/ NOW | Nw" |^5cMpvZ\A 


AMT_ DIR. 

ARA 

CCP-GAIN 

CODE 

CON-DELTAV 
DELAY 
DO DC 
FILTER 
G£MM 
BOW- FACTOR 


PRP ! 

PW-BAUD 1 %ms 

RX_ DELAY 
RX-GAIN 
SAMPLE- MIN 
SPACE I ! U 

SPAN 
SPECTRA 
SSTDODC 
TDA 


WINDOW 

DOPOWCOR. 


I ml~nki\ 





















ETEE5 


CONDO: 




DATE OF CHANGE: Jgj'f*?, csr 


PAEAMCTCK 


ANT-DIE 

AEA 

CCP-GAIN 

CODE 

'CON-DELTAY 


DELAY 

DODC 

FILTER 

GCNUM 

FOW-FAETOR. 


PRP 

PW.BAUD 
RX_ DELAY 
EX-GAIN 
SAtoPLE-MlN 


SPACE 

SPAN 

SPECTRA 

SSTDODC 

TDA 


WINDOW 

DOPOWCOR 


NODE 1 


UWIT I NORM | NEW 



Dt 2> KAQDE4 I N0DC 5 1 ^ODE16 


mdbmI new I noumnon || rcsm I new I noen| new 


COMMON PAR MS |UNU 


BLANK-TX 

coM5£sus.ser 


COOLDOWN 

FREED6LAY 

NUNV6ATE5 

NUNlN^CfcS 

•TR.TX 


TX- BLANK. 

TC-DUTY 

TX-TR 



[ 

S s -«4* 

!i 

MU 4 

T 

fax? 

y 































>AR PARMIEES 

11 r 

rmt-ncj-. gw (.c^DKlOF CHW6C: WiUcv/??/ H ■ «<> 


PA.ZAMCTE& 


MODE \ 


UMIT I MOAtt I NEW 


MODE 3 I MODE 4 II MODE 5 I MODE 6 


NOAM I NEW 1 ADA* I NEW 


Mew II PCKIA I NEW I NOV* I NEW 


AKJT_D1£ , 

ara 

CC.P-GAIN 

CODE 

CON.DQXAV 


DELAY 

DODC 

FILTER 

GCNUM 

ROW-FACTOR. 


PRP 

PW. BAUD 
RX- DELAY 
EX -GAIN | /V 
SAMPLE- NUN /V 


SPACE 

SPAN 

SPECTRA 

SSTDODC 

TDA 


WINDOW 

DOPOWCOR 


I i wRga I 


rfVl 




WWi 


fVJMOAl PARMS lawn 


BLANK-TX 

CON5E5US-SET 


COOLDOWN 

FREEDELAY 

NUMGATES 

NUMMODES 

TR-TX 


TX-BLANK 

TY-DUTY 

TX-TR 


( > H^o -e 

H~o A/ e S «>Hr* > 

































l\f\0C£7 I (\A0D6 3 I N\DDt 9 1 MOO. <0 | W006 H [MOE \Z 


PARAMETER II UNIT II UWrtfl NEW I NORM 


NDEMB NEW R NOSH I NOW H NOW* |NE\W 



WV WZTu ~ , ^ q 'i *7 r’H' 

























"i"iw ifpiMiwi mi' m 




